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Abstract 

The proton-nucleon cross section ratio R = Br(T — • da{T) / dy\y=Q / a{J/ijj) 
has been measured with the HERA-B spectrometer in fixed-target proton-nucleus 
collisions at 920 GeV proton beam energy corresponding to a proton-nucleon cms 
energy of y/s = 41.6 GeV. The combined results for the decay channels T — > e'''e~ 
and T jJ''^ IJ.~ yield a ratio R = (9.0ib2.1) • 10~^. The corresponding T production 
cross section per nucleon at mid-rapidity [y = 0) has been determined to be Br(T — ^ 
/+r ) • da{l:)/dy\y=Q = 4.5 ± 1.1 pb/nucleon. 

Key words: 

T mesons, cross section, proton-nucleus collisions, -y/s = 41.6 GeV 
PACS: 13.20.Gd, 13.85.Qk, 14.40.Gx, 24.85.+p 
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1 Introduction 



In recent years, there has been a rapid development of models describing 
quarkonium, especially charmonium, production, with great success both in 
the high- and low-energy range (see e.g. review [1,2] and references therein). 
These developments are driven by the available measurements of charmonium 
production, and the implication of a possible suppression of charmonium pro- 
duction as an indicator for a Quark-Gluon Plasma (QGP) [3,4]. Bottomonium 
production represents a natural field for testing the predictions of these mod- 
els [1] . Measurements of the T production cross section ^ have been performed 
by many experiments [5]- [18] in the wide range of the proton-nucleon centre- 
of-mass (cms) energy a/s of 19 to 1800 GeV, and with targets ranging from 
proton {A — 1) to platinum {A — 195) with both proton and antiproton 
beams. 

Usually, the result of bottomonium production measurements is presented as 
the product of the differential cross section at mid-rapidity da{T)/dy\y^Q times 
branching ratio [1]. However, the acceptance of most of the fixed-target ex- 
periments which have so far measured bottomonium production, is not in the 
region of central collisions, y = 0, so that systematic effects in determining the 
total cross section can be substantial. More specifically, precisely in the en- 
ergy region of HERA-B, the available results [10,12] disagree by about a factor 
of two (see Table 1 ^ ) which cannot be explained by a large nuclear suppres- 
sion [19]. HERA-B covers the region of mid-rapidity, implying less uncertainty 
in the determination of the total cross section, and can thus contribute to the 
knowledge of T production, despite its rather small sample size. Furthermore, 
both T — > and T — > e+e~ decay channels are measured simultaneously, 

providing an additional statistically independent cross check. 



2 Measurement method 

Wc determine the T production cross section by comparing the relative yields 
of T and J/ ip production, and normalizing to the known J/ ip cross section: 

^ Throughout this paper we refer to the sum of T(IS') -|- T{2S) + T(3S') cross 
sections as the T cross section. 

^ Both papers [10,12] quote the differential cross section in terms of the Feynman 
scaling variable xp which we transform to rapidity y. The two quantities are related 
via d(T{T) / dxp\^^^Q = F{^/s) ■ da{T)/dy\y=o, where F{^/s) is a coefficient which 
depends on the (measured) transverse momentum distribution of the produced T 
mesons. Its numerical value is 1.98 ± 0.03 and 2.12 ± 0.03 for ^/s = 38.8 GeV and 
■s/s = 41.6 GeV, respectively. 
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Table 1 

Summary of the available measurements of the T production cross section in pA 
collisions near ^/s = 41.6 GeV. The published result of [12] refers to the production 
of T(1S) only and has been rescaled using Eq. (4) to also include T(2S) and T(3S). 
In the case of E771 [13], the published value has been corrected to include the T(3S) 
state and the coefficient At/efr (see Sect. 4.4). Overall normalization uncertainties of 
15% (Ref. [10]) and 10% (Ref. [12]) have been included. 



GeV 


BTiT^l+l-).da{T)/dy\y^„ 
pb/nucleon 


Tar- 
get 


Expe- 
riment 


38.8 


2.11 ±0.33 


Cu 


E605 [10] 


38.8 


2.31 ±0.38 


Be 


E605 [11] 


38.8 


4.7 ±0.5 


D 


E772 [12] 


38.8 


7.7 ±3.2 


Si 


E771 [13] 



Br(T^Z+r)-^(T) 
dy 



and in addition, we define the ratio of the mid-rapidity T cross section to the 
total J/ip cross section as 



Br(T ^ 1+1-) • da{T)/dy\y=o 
a{Jli^) 



(2) 



Here, a{J/ip) is the J/ip production cross section, and 7V(T) and N{J/ip) are 
the numbers of observed T and J/ if) decays, respectively. The branching ratio 
J/V^ ^ is taken as the average of the most recent values for 7/^0 —> e^e~ 
and J/ip ^ fi^fi^ [20], e{J / ip) / e{T) is the ratio of J/ip and T trigger and re- 
construction efficiencies determined from Monte Carlo simulations, and Ayefr is 
the coefficient relating the full and differential T cross section at mid-rapidity 
y — {): (j(T) = A^/eff da{T) / dy\y^Q (see Sect. 4.4). In this way the result is in- 
dependent of the luminosity determination, and systematic uncertainties due 
to absolute trigger efficiencies cancel out to a large extent, since only relative 
efficiencies between T and J /ip and relative acceptances enter the final result. 
For a{J/ip) at ^/s = 41.6 GeV we use the value 



<^pn{J/iP) = 502 ± 44 



nb 



nucleon 



(3) 
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obtained from a global fit [21] of J /if} production data^, after adjusting the 
data using a common nuclear suppression parameter a = 0.96 ±0.01 [19] and 
the latest value of branching ratio Br^J/ip — > [20]. Eq. (1) assumes the 

same numerical values of a parameters for T and J/ip production in accor- 
dance with the experimental data [19,23]. 



3 The HERA-B detector and the data sample 

HERA-B is a fixed-target experiment at the 920-GeV HERA proton storage 
ring of DESY and consists of a forward magnetic spectrometer featuring a 
high resolution vertexing and tracking system and offering a good coverage 
of the central region of collisions (the Xp range is about [—0.6,0.15] for T 
production and about [—0.35,0.15] for J/ip production, corresponding to the 
rapidity ranges about [—1.2,0.3] and [—1.5,0.5], respectively). 

The main components of the detector are sketched in Fig. 1. The target con- 
sists of wires of various materials which are inserted into the halo of the HERA 
proton beam. Data are taken with carbon, tungsten and titanium target wires 
operated at an interaction rate between 5 and 8 MHz. The Vertex Detec- 
tor System (VDS) consists of silicon micro-strip detectors located within the 
vacuum vessel of the target. The first station of the main tracker is placed 
upstream of the 2.13 T-m dipole magnet and the remaining 6 tracking sta- 
tions are placed downstream. Muon identification is performed by the muon 
detector (MUON) , while the electrons are detected and identified by the elec- 
tromagnetic calorimeter (ECAL). 

The trigger chain includes pretriggers provided by ECAL and MUON for the 
lepton candidate search, and a first level trigger (FLT) which finds tracks 
downstream of the magnet starting from the pretrigger seeds. The FLT re- 
quires that at least two pretrigger candidates be present in an event and that 
an FLT track be found from at least one of them. The (software) second level 
trigger (SLT) starts from pretrigger candidate tracks, confirms them in the 
tracker and VDS using a simplified Kalman filter algorithm, and accepts the 
event if either two electron or two muon candidates with a common vertex are 
found. The trigger imposes a cut on the transverse energy Et of the electrons, 
and an implicit cut on the transverse momentum pt of the muons. A more 
detailed description of the HERA-B detector, trigger and reconstruction chain 
can be found in Refs. [24,25], and references therein. This analysis is based on 

^ The global fit includes our own measurement of the J/ip cross section 663 ib 74 ± 
46nb/nucleon [22]. We prefer to normalize the present measurement to the value 
given by the global fit since the fit provides a complete summary of all available 
measurements. 
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Om 5 10 15 20 



Fig. 1. Plan view of the HERA-B detector. 
134 ■ 10^ events obtained in the 2002-2003 physics run. 



4 Data analysis 

4-1 Monte Carlo simulation, trigger and reconstruction ejficiency 

Monte Carlo simulation is performed for J/ip, T and Drell-Yan production. 
Heavy quarkonia and Drell-Yan production are simulated with PYTHIA 5.7 
[26]. The energy remaining after the quarkonium or Drell-Yan generation is 
passed to FRITIOF 7.02 [27] to simulate the underlying pA interaction. Fi- 
nally, additional inelastic pA interactions, generated by FRITIOF, are super- 
imposed according to the multiple-interaction probabilities of five separate 
running periods. After tracking the particles through the detector material 
with GEANT 3.21 [28] and a realistic digitization simulation, the event is 
reconstructed by the same reconstruction program which is used for the real 
events. The trigger and the detector parameters are tuned for the individual 
data taking periods. To have a realistic description of the kinematic distribu- 
tions for J/%1) and T, the PYTHIA generation has been tuned by reweight- 
ing according to the differential distributions from the high-statistics data of 
Refs. [10,12,29,30]. The relative trigger and reconstruction efficiencies are de- 
termined from Monte Carlo simulations. For the muon channel, the ratio of 
efficiencies is e{J / ip) / e{T) = 0.76 ±0.05. The departure from unity is a result 
of the trigger px cut which affects muons from J/ ip decay more strongly than 
muons from T decay due to the different mass scale. Because of a harder cut 
on Et at the trigger level, the difference is even larger in the electron channel: 
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£(J/V')/£(T) = 0.31 ±0.03. 



4-2 Event selection 

The event selection includes cuts on the quality of tracks, most notably, that 
the tracks must have segments in both the VDS and the main tracker and 
that "clones" (nearby reconstructed tracks originating from the same real 
physical track) be removed. Further requirements are that either two muon 
or two electron candidates of opposite charge and with a common vertex are 
present. Muon candidates are required to penetrate to the muon detector layers 
behind the absorber material. Electrons are identified by requiring an energy 
deposit (cluster) in the ECAL. Bremsstrahlung photons emitted by electrons 
in front of the magnet are recovered for better energy resolution and additional 
electron identification [25]. The transverse energy Et at ECAL must exceed 
0.95 GeV, and the ratio between the energy as measured by the ECAL and 
the momentum must be close to unity (0.85 < E/p < 1.25 and, in addition, 
E/p > 0.9 when no bremsstrahlung photon is found). 

To suppress background from secondaries produced in the beam pipe which 
reach the muon chambers, we impose cuts on the muon transverse momentum, 
Pt, which are linear functions of the mass with cut ranges: pt G [3, 6] GeV/c 
for the invariant mass m — mr, Pt £ [0.7,4] GeV/c for m — rrij/^ and 
cuts on the muon total momentum p: p & [15,190] GeV/c for m = mx, 
p G [4,100] GeV/c for m = rrij/^; for intermediate dilepton masses a lin- 
ear interpolation of the cut limits is used. For electrons which do not suf- 
fer from a similar background, the cut for the electron total momentum is 
p e [4,190] GeV/c and for the transverse momentum pr G [0.7,6] GeV/c. 
After applying these cuts, the Icpton momentum spectra for the mixture of 
signal and all background sources are similar in the data and Monte Carlo. 

4-3 Description of signal and background 

4-3.1 Muon channel 

The dimuon mass spectrum in the J/ ip mass region (Fig. 2) is fitted with a 
Gaussian folded with contributions from J/ip ^ /x+/x~7 [31] for the signal 
shape, and an exponential of a second-order polynomial for the background. 
For the combined sample of the three target materials (64% C, 33% W, 3% Ti) 
we obtain a total of 152900±490 J/ip decays and a width of 38.8±0.1 MeV/c^. 

For the mass region m > 5 GeV/c^ (Fig- 3) we use a similar function to describe 
the T peaks to that used for the description of the J/t/j. Due to the lack of 
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r I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

2.6 2.7 2.8 2.9 3 3.1 3.2 3.3 3.4 3.5 

11^ ji invariant mass, GeV/c^ 

Fig. 2. Fit of the dimuon mass spectrum obtained after J/V^ selection cuts. 



statistics, the positions of the T{2S) and T(3S') states are fixed relative to the 
T(IS') state using the PDG mass values [20], and the relative contributions of 
the three states are fixed according to the E605 results [10]: 

N{IS) : N(2S) : N(3S) = (70 ± 3) : (20 ± 2) : (10 ± 1). (4) 



The width of the T(IS') state is fixed to the width of the J/ip scaled by the 
ratio of the expected momentum resolution for muons from T and J/ ip decays 
resulting in 159 MeV/c^. The widths of the T(25') and T(35') states are scaled 
proportionally to their masses. 

In this mass region, the Drell-Yan process as well as random combinations be- 
tween leptons contribute to the background. As before, we describe the com- 
binatorial background by an exponential of a second-order polynomial, whose 
shape is determined from a fit to the like-sign yU^/i^ pair spectrum (Fig. 4), 
and a normalization factor which is left as a free parameter. The shape of 
the Drell-Yan spectrum is determined from a fit to the corresponding recon- 
structed Monte Carlo events. 

Thus, the free fit parameters for the spectrum with m > 5 GeV/c^ (Fig. 3) 
are: the total T yield, the T(15') mass, the yield of Drell-Yan dileptons, and 
the height of the combinatorial background. The total T yield in the log 
likelihood fit is N{T) = 30.8 ± 7.4stat- The fit function describes the data well 
(xVndf = 51/54). Using this value in Eq. (1) we obtain R'^'J^ = (7.9 ± 1.9stat) • 
10^®, and using the value of the J/ip cross section Eq. (3), the cross section 
for T production at mid-rapidity as determined from T — > iJ,~^n~ becomes 
Br(T ■ da{T) / dy\y=o — 4.0 ± l.Ogtat pb/nucleon. 
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Fig. 3. Fit of the T signal (thick line) with the individual contributions of 

the T(15'), T{2S), T(3S) states shown in thin solid lines. The background consists 
of Drell-Yan pairs (dashed line) and combinatorial contribution estimated from the 
H^H^ spectrum (dotted line). 




5 6 7 8 9 10 11 12 

/i^jU^ invariant mass, GeV/c^ 

Fig. 4. Spectrum of like-sign /x^/x^ pairs, together with a fit to an exponential of a 
second-order polynomial. 

The fit (Fig. 3) shows that the combinatorial background dominates below 
7 GeV/c^ but becomes negligible at the T mass. We have verified that the 
generated shape of the Drell-Yan background in Monte Carlo simulations 
agrees with the published data [10,12,32]. Further support for the correctness 
of the background curves in Fig. 3 comes from studying a number of sensitive 
kinematic variables (such as the average momentum asymmetry between the 
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muons). The observed values are compatible with the mixture of signal and 
background contributions obtained from the fit to the mass spectrum. 



The fit also provides a measurement of the Drell-Yan (DY) differential cross 
section (Pa /{dmdy)\y=Q in the region of the T mass, or equivalently the ratio 
i?DY [9,10,11] of the T production to the Drell-Yan cross sections 



da 



i^.. = Br(T^/-r).-(T)^^J^(DY) 



y = 
m=9.46 



Similarly to Eq. (1) we determine for the Drell-Yan cross section: 
£^y{^^)\llL = 2.0 ± 0.5pb/(micleon.GeV/c2) and R^^ = 2.0 ± 
0.8 GeV/ where statistical and systematic uncertainties have been combined. 

To compare these results with the published data of the Fermilab fixed-target 
experiments [10,11,12,32] at ^/s = 38.8 GeV, we scale the Drell-Yan pro- 
duction cross section with the variable r = m^/s leading to a factor of 
1.28 ± 0.03^, and the T cross section by a factor 1.32 (obtained from the 
fit to the T data, see Sect. 5). The results, scaled to ^/s — 38.8 GeV, then be- 
come: £^{DY)\f£o^^^ = 1.5 ± 0.4pb/(nucleon • GeV/c^) and i?^^|^=38-8 = 
2.0 ± 0.8 GeV/c^, which agrees within the uncertainties with the fits of the 
data from the Fermilab experiments: i?DY between 1.3 and 1.6 GeV/c^ [10,11] 
and ££-(DY)|^:j;,g = 1-4 ± 0.3pb/(nucleon • GeV/c^) [10,12,32]. 



4-3.2 Electron channel 

The analysis of dielectron events generally follows the path described for 
dimuon events. Before fitting, the momentum vectors of the leptons are cor- 
rected by adding in the energy of bremsstrahlung photons emitted in the 
material before the magnet and reconstructed in the calorimeter. The fit 
function takes into account resolution effects as well as a tail due to non- 
recovered bremsstrahlung and final state radiation [25,31]. The fit yields 
NiJ/ijj) = 109710 ± 930 and a width of 59.9 ± 0.8 MeV/c^ for the combined 
sample of the three target materials (62% C, 32% W, and 6% Ti). 

The log likelihood fit of the mass region m > 5 GeV/c^ (Fig. 6), for which 
the free parameters are the position of the T(IS') peak, the contributions of 
the T signal, of Drell-Yan, and of the combinatorial background, results in 
75 ± 14stat events with x^/ ndf = 19/27. This leads to a cross section ratio 
Rj^/^ — (11.0±2.1stat) • 10~^, and a mid-rapidity cross section for T production 
of Br(T — > e+e") • da{T) / dy\y^Q — 5.5 ± l.Ogtat pb/nucleon. 

^ This number agrees well with the theoretical prediction [33] of about 1.275 for 
the mass region 8 < m < 11 GeV/c^ . 
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Fig. 5. Fit of the dielectron mass spectrum obtained after J/ if; selection cuts. 




e+e invariant mass, GeV/c^ 



Fig. 6. Fit of the T — e+e" signal (thick line) with the individual contributions of 

the T(15), T(2S'), T(35) states shown in thin solid lines. The background consists 
of Drell- Yan (dashed line) and combinatorial contributions (dotted line) . 

As for the muon channel, the various energy-dependent kinematic distribu- 
tions of the electron pairs behave as expected from the relative contributions 
of the different processes. The imposition of an Et requirement for electrons 
in the trigger severely suppresses triggers from like-sign leptons, because of the 
additional /^^-bias of the magnet. The resulting like-sign electron mass spec- 
trum therefore does not describe the combinatorial background and is not used 
in this analysis. Instead, the combinatorial background is evaluated by mix- 
ing tracks from different events. For masses greater than about 6 GeV/c^, the 



12 



combinatorial background is less than the Drell-Yan contribution and becomes 
neghgible for masses above 9 GeV/c^. 



The Drell-Yan contribution at m = 9.46 GeV/c^ is 



cPa 
dmdy 



(DY) 



y = 
m=9.46 



3.7 ± 1.0pb/(nucleon • GeV/c^) leading to R^y = 1.5 ± 0.7 GeV/c^. Scal- 
ing these results to ^/s — 38.8 GeV we obtain 



y = ^ = 2.9 ± 

m=9.46 



dm dy ("^"^) 

0.8pb/(nucleon-GeV/c2), and R'^-y\^'=^^-^ = 1.4 ± 0.6 GeV/c^. The Drell 
Yan cross section measurements for the electron and muon cases differ by a 
factor of 1.9, corresponding to 1.6 standard deviations; the ratio R^y a-grees 
with the muon case and with Fermilab measurements. 



4-4 Systematic uncertainties 



The systematic uncertainties are dominated by the uncertainties in the de- 
scription of the background which contribute 14% for muons and 17% (in- 
cluding uncertainties in the bremsstrahlung tail) for electrons, respectively. 
The uncertainties of the relative J/t/j and T efficiencies are estimated to be 
7% for muons and 9% for electrons. The systematic uncertainty of the J/ijj 
reference cross section Eq. (1) is 9%. The parameter Aycs{^/s) = 1.14±0.12syst 
(precision of 11%) at ^/s = 41.6 GeV is determined from the fits of existing 
measurements for xp and px distributions for T mesons [10,12]. 

Other systematic uncertainties are small compared with those already men- 
tioned. Among them arc the use of a Gaussian function with a fixed width 
for fitting of the T peaks (< 4%), uncertainties of the fractions of the various 
T states in Eq.(4) (< 0.1%), the polarization effects in J/t/j and T produc- 
tion (< 1.8%), and the branching ratio J/ip l^l" (< 1.7%). When, in 
the electron channel, the particle identification requirements are strengthened 
by requiring either that one of the two lepton candidates have an associated 
bremsstrahlung cluster in the calorimeter, or that both have associated clus- 
ters, the results change by less than 4%. The results are stable within the 
statistical uncertainty for a wide variation of the cuts for muon or electron 
identification. All contributions, added in quadrature, result in a systematic 
uncertainty of 21% in the case of muons, and 25% in the case of electrons. 



5 Combined results 



Table 2 summarizes all results obtained in the previous sections. One can see 
that the results in the muon and electron channels are compatible. Combining 
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Table 2 

Summary of the input vahics and the res 
clccLrou channels. 



cross sections for both muon and 



Parameter 


fi'^fj, channel 


e^e channel 


N(T) 


30.8 ± 7.4 


75± 14 


NUIxh) 


152900 ± 490 


109710 ± 930 


£(JM)/£(T) 


0.76 ± 0.05 


0.31 ±0.03 


(ypN{J/i') 

(nb/nucleon) 


502 ± 44 


Br(J/V' ^ Z+r) 


(5.90 ± 0.10)% 


Ayeff 


1.14±0.12,yst 


Rjl^ (units of 10-6) 


7.9±1.9stat±1.5syst 


11.0 ± 2.1stat ± 2.5syst 


Br(T-^/+/-)-^(T)l . 
^ dy '^y-^ 
(pb/nucleon) 


4.0 ± l.Ostat ± O.Ssyst 


5.5±1.0stat±l-4syst 


dV(DY) 

dmdy 
(pb/nucleon • GeV/c^) 


2.0 ±0.5 


3.7 ±1.0 


Rdy 


2.0 ±0.8 


1.5 ±0.7 



the muon and electron channels, we have 



R 



(9.0 ±2.1) • 10" 



and, using the J/t/j reference cross section Eq. (3), we obtain (Fig. 7) 



Br(T^/+r)-^(T) 
dy 



= 4.5±1.1^^^ 
2/=o nucleon 



(5) 



where the error includes both statistical and systematic contributions. The 
systematic uncertainties in the muon and electron channel have a common 

part of 14% compared to the full values of 21% and 25%, respectively, the 
dominant contributions being due to the uncertainties of Aycs and the J/ip 
reference cross section Eq. (3). After taking this into account, the of the 
combined result is — 0-6. 

If, instead of normalizing to the J/ip cross section, we use the yield 
and efficiency ratio between T and Drell-Yan production together with 
the DY cross section at V« = 38-8 GeV of £i^(Dy)\y=o = (1-35 ± 
0.27) pb/(nucleon • GeV/c^) [10], scaled by the factor 1.28 ± 0.03 to 
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20 25 30 35 40 45 50 55 60 65 

Vs, GeV 

Fig. 7. This result (HERA-B) compared to the fit (solid curve) of the data of pre- 
vious experiments [5]-[16] up to = 63 GeV. The NLO theoretical curves (see 
description in the text), shown as a dashed line, differ by less than the width of the 
line. 



^/s = 41.6 GeV, we obtain for the T production cross section: 
Br(T ^ m-) ■ ^(T) = 2.9 ± 1.1 with a of 0.3. However, the 

^ ' ay ^ ' y—Q nucleon ^ ' 

pubhshed data on the DY process are less abundant than those for the J/^ 
cross section and are less precise due to limited statistics and difficulties in 
separating the T signal from the DY continuum. 

We therefore take the value obtained by normalization to the DY cross section 
as a confirmation of our main result Eq. (5), which is obtained by normalization 
to the J/t/j cross section. Scaling this result by a factor of 1.32 (obtained from 
the fit described below) to compare it to the published data at ^/s = 38.8 GeV, 
we obtain: Br(T ^ /+/") • g(T)|^=^^-^ = 3.4 ± 0.8 pb/nucleon. Our result 
lies half-way between the results of E605 [10,11] and E772 [12] and does not 
favour one of these results over the other. In Fig. 7 we compare the result 
Eq. (5) with the measurements of previous experiments. In most cases, the 
uncertainties are large with the exception of four measurements obtained at 
Fermilab at cms energy ^/s = 38.8 GeV, which are also those closest in energy 
to HERA-B {y/s = 41.6 GeV). However, these values are in poor agreement 
with each other, as summarized in Table 1. 

The fit of all experiments in Fig. 7 (solid line) uses the Craigie parameterization 
[34] /(^/s) = (Jo exp (— mo/v^), yielding ctq = 182 ± 21 pb/nucleon, uiq = 
167±4 GeV, x^/ndf = 37/14. The dotted line in Fig. 7 shows fits to predictions 
of next-to-leading order (NLO) calculations from Ref. [35] in the framework of 
the colour evaporation model (GEM) using the MRST HO parton distribution 
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functions [36]. The thickness of the hne corresponds to the variations of results 
in the three sets of input parameters of the model which describe the open 
beauty production data [35]: nij, = fi = 4.75 GeV/c^; rrib = 4.5 GeV/c^, 
/i = 2mf,; and rrih = 5 GeV/c^, /i = mh/2, where is the 6-quark mass and fj, 
is the renormalization scale (assumed to be equal to the factorization scale). 
The NLO predictions are normalized according to a fit of the experimental 
data using the mentioned values of the parameters and fi [35]. 



6 Conclusion 



The T production yield at mid-rapidity in pA collisions at a proton momentum 
p = 920 GeV/c has been measured in both channels T — ii~ and T — *■ 
e^e^ . The J / ip cross section (3) has been used for normalization, and the ratio 
of T and J /ip cross sections Eq. (2) is determined to be 

Rj/^ = (9.0 ± 2.1) ■ 10-^ 

The resulting T production cross section (both lepton channels combined) is 



Br(T ^ Z+r) ■ ^(T) 
dy 



4.5 ± 1.1 



pb 



y=o nucleon 



Our result, scaled for ^/s dependence, lies half-way between those of E605 
[10,11] and E772 [12]. Normalization with respect to the Drell-Yan process 
confirms this result. The result agrees within 1.4 standard deviations with 
current GEM NLO predictions [35] (see Fig. 7). 



Acknowledgments: 

We express our gratitude to the DESY laboratory and to the DESY accel- 
erator group for their strong support since the conception of the HERA-B 
experiment. The HERA-B experiment would not have been possible without 
the enormous effort and commitment of our technical and administrative staff. 
It is a pleasure to thank all of them. We would like to thank R. Vogt for many 
stimulating discussions and suggestions. 



References 



[1] R. Vogt, Phys. Kept. 310 (1999) 197. 

[2] N. Brambilla et al., CERN Yellow Report CERN-2005-005, |hep-ph/0412T58l 



16 



[3] J.F. Gunion and R. Vogt, Nucl. Phys. B492 (1997) 301. 

[4] F. Karsch, M.T. Mehr and H. Satz, Z. Phys. C37 (1988) 617. 

[5] NA3 collab., J. Badier et al., Phys. Lett. B86 (1979) 98. 

[6] J.K. Yoh et ah, Phys. Rev. Lett. 41 (1978) 684; Erratum ibid. 41 (1978) 1083. 

[7] W.R. Innes et al., Phys. Rev. Lett. 39 (1977) 1240; Erratum ibid. 39 (1977) 
1640. 

[8] K. Ueno et al., Phys. Rev. Lett. 42 (1979) 486. 

[9] S. Childress et al., Phys. Rev. Lett. 55 (1985) 1962. 



[lo; 

[11 
[12 

[13; 

[14 
[15 

[16; 
[ir 

[is; 

[19 
[20 
[21 
[22 

[23; 
[24 

[25 

[26; 

[27 

[28; 



E605 collab., G. Moreno et al., Phys. Rev. D43 (1991) 2815. 

E605 collab., T. Yoshida et al., Phys. Rev. D39 (1989) 3516. 

E772 collab., P.L. McGaughey et al., Phys. Rev. D50 (1994) 3038; Erratum 
ibid. D60 (1999) 119903. 

E771 collab., T. Alexopoulos et al., Phys. Lett. B374 (1996) 271. 

R-108 (CCOR collab.), L. Camilleri, in: T.B.W. Kirk, H.D.I. Abarbanel (Eds.), 
Proc. 1979 Int. Symp. on Lepton and Photon Interactions at High Energies, 
Fermilab (1979), p. 282. 

C. Kourkoumelis et al., Phys. Lett. B91 (1980) 481. 

CCOR collab., A.L.S. Angelis et al., Phys. Lett. B87 (1979) 398. 

UAl collab., C. Albajar et al., Phys. Lett. B186 (1987) 237. 

CDF collab., F. Abe et al., Phys. Rev. Lett. 75 (1995) 4358. 

E866/NuSea collab., M.J. Leitch et al., Phys. Rev. Lett. 84 (2000) 3256. 

S. Eidelman et al., Phys. Lett. B592 (2004) 1. 

F. Maltoni et al., |hep-ph/0601203| (2003), submitted to Phys. Lett. B. 

HERA-B cohab., I. Abt et al.. Preprint DESY 05-232, Hamburg (2005), 
,hep-ex/0 512029, submitted to Phys. Lett. B. 

E772 collab., D.M. Aide et al., Phys. Rev. Lett. 64 (1990) 2479. 

HERA-B collab., E. Hartouni et al., HERA-B Design Report, DESY-PRC 95-01 
(1995); HERA-B collab., HERA-B Status Report, DESY-PRC 00-04 (2000). 

HERA-B cohab., I. Abt et al.. Preprint DESY 05-233, Hamburg (2005), 
|hep-ex/0512030^ Phys. Rev. D (in press). 

T. Sjostrand, Comp. Phys. Comm. 82 (1994) 74. 

H. Pi, Comp. Phys. Comm. 71 (1992) 173. 

R. Brun et al., GEANT3, CERN-DD-EE-84-1 (1987). 



17 



[29] E789 collab., M.H. Schub et al., Phys. Rev. D52 (1995) 1307; Erratum ibid. 
D53 (1996) 570. 

[30] E771 collab., T. Alexopoulos et al., Phys. Rev. D55 (1997) 3927. 

[31] A. Spiridonov, Preprint DESY 04-105, Hamburg (2004), plep^ex/0510076[ 

[32] E866/NuSea collab., J.C. Webb et al., pIep^x/0302019 (2003); 
J.C. Webb, hep-ex/0301031 (2003); 

E866/NuSea Collab., C.A. Gagliardi et al., AIP Conf. Proc. 698 (2004) 100. 
[33] R. Vogt, private communication (Feb. 2005). 
[34] N.S. Craigie, Phys. Rept. 47 (1978) 1. 
[35] M. Bedjidian et al., |hep-ph/031 10 481 (2003). 
[36] A.D. Martin et al., Eur. Phys. J. C28 (2003) 455. 



18 



